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SBIR  Final  Report  -  GIS  for  Marine  Operations 
Executive  Summary 

An  SBIR  Phase  I  project  was  undertaken  by  Delta  Data  Systems,  Inc.  and  Neptune 
Sciences,  Inc.  which  focused  on  feasibility  issues  relating  to  the  development  of  a 
Marine  Geographic  Information  System  (MGIS)  that  could  support  Logistics  Over  the 
Shore  (LOTS)  planning,  operational  monitoring  and  analysis  tasks.  Significant 
feasibility  issues  relating  to  MGIS-numerical  model  interaction  were  addressed  in  the 
Phase  I  effort. 

Our  view  of  MGIS  support  to  LOTS  involves  visualization  of  the  coastal  environment 
and  coastal  processes  from  both  an  historical  and  a  contemporary  (near  real-time) 
perspective.  MGIS  can  organize  information  in  a  geospatial  context.  Making  the 
context  sufficiently  "information  rich"  to  support  interactive  decision  making  is  a  function 
of  data  availability  and  quality. 

We  believe  that  GIS  is  a  mature  technology  that  has  established  an  ability  to  create, 
integrate,  display  and  analyze  georeferenced  data  in  a  visually  informative  system. 
However,  applying  GIS  visualization  and  analysis  to  planning  problems  in  the  littoral 
environment  poses  challenges  in  data  collection,  data  integration,  visualization  and 
mapping  that  are  generally  not  encountered  in  land  based  GIS  projects. 

Unlike  land  based  GIS  development  which  (in  several  industrialized  countries)  is 
supported  by  ready  access  to  a  well  developed,  relatively  high  resolution,  spatial  data 
infrastructure,  marine  based  GIS  development  must  rely  on  comparatively  under¬ 
developed,  low  resolution  geophysical  and  meteorological  data  bases.  In  addition,  the 
dynamic  nature  of  the  near  shore  and  surf  zone  "landscape"  requires  that  much  of  our 
view  must  be  developed  through  modeling.  Modeling  relies  on  hindcasts,  forecasts 
and/or  near  real  time  observations  of  sea  state  and  weather  conditions  (in  combination 
with  physiographic  data)  to  develop  an  estimate  of  their  combined  effect  on  (normally  a 
very  small)  surface  region  near  the  coast  including  the  surf  zone. 

The  principal  challenges  in  determining  the  feasibility  of  and  planning  the 
implementation  of  MGIS  were: 

1 )  Data  collection  and  data  integration. 

Efficient  and  effective  strategies  must  be  designed  to  enable  rapid  location  and 
integration  of  geophysical  and  meteorological  data  of  two  general  types,  historical  and 
near  real-time  with  different  characteristics. 


A  survey  of  data  sources  and  marine  information  systems  was  completed  in  Phase  I. 
Navoceano's  IDBMS,  the  Master  Environmental  Library  and  the  Satellite  Active  Archive 
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were  identified  as  examples  of  "area  locator"  information  systems  which  represent  new 
wave,  internet-based  tools.  In  addition,  geoprocessing  tools  were  applied  to  the 
development  of  high  resolution  surf  zone  bathymetry  and  shoreline  as  a  demonstration 
of  critical  data  development  in  cases  where  no  digital  data  are  available. 

2)  Integration  of  numerical  modeling  with  GIS. 

A  workable  exchange  between  selected  numerical  models  and  MGIS  must  be 
developed.  The  specific  concerns  are  how  to  use  geoprocessing  tools  to  "feed" 
environmental  and  geophysical  data  to  models  and  how  to  integrate  and  visualize 
model  output  in  a  geographic  context  for  application  to  LOTS  planning  and  operations. 

As  an  example.  Phase  I  successfully  demonstrated  an  interchange  between  the  Navy 
surf  model  and  a  conceptual  MGIS.  The  surf  model  was  "fed"  by  MGIS  and  surf  model 
output  was  visualized  as  a  map  of  predicted  wave  and  surf  conditions. 

3)  Human  Factors. 

A  task  oriented,  graphical  user  interface  must  be  created  that  effectively  supports 
varying  levels  of  user  background  and  expertise. 

The  development  of  a  conceptual  MGIS  in  Phase  I  provided  insight  into  task 
identification.  Basic  software  design  concepts  for  a  Phase  II  implementation  were 
developed.  A  concept  MGIS  was  developed  and  delivered  to  the  SBIR  Technical 
monitor  at  WES  that  included  data  and  representative  geoprocessing  tools. 
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Project  Overview:  Objectives  and  Results 
Objectives: 

1 .  To  identify  technical  feasibility  issues  relating  to  the  integration  of  geoprocessing 
and  numerical  modeling  to  support  LOTS  planning  and  operations.  Specific  focus 
was  on  MGIS-numerical  model  interaction.  To  demonstrate  an  exchange  of  data 
and  information  between  a  geoprocessor  and  a  numerical  model,  we  developed  a 
test  which  would  generate  and  then  visualize  a  surf  zone  forecast  in  a  detailed 
geographic  context.  (Although  several  models  were  identified,  the  Navy  SURF 
model  was  the  only  surf  forecasting  software  available  for  operation  during  Phase  I.) 

2.  To  investigate  sources  of  geophysical,  oceanographic  and  meteorological  data  in 
digital  form  and  to  assess  their  accessibility  and  utility  to  a  LOTS  planning  and 
decision  support  package.  Again,  the  primary  focus  was  on  data  related  to  the 
requirements  of  a  surf  forecasting  model.  Foremost  among  these  were  bathymetry 
and  observations  of  wave,  wind  and  tide.  This  part  of  the  effort  was  performed 
specifically  to  provide  insights  into  LOTS  planning  and  decision  support  design 
issues. 

3.  To  test  various  "processing  hypotheses"  relating  to  data  generation  and 
visualization.  Specifically,  we  were  interested  in  the  applicability  and  adaptability  of 
an  existing  set  of  surface  modeling  and  profiling  tools  to  the  generation  of 
bathymetry  and  the  extraction  of  numerical  model  input  parameters. 

4.  To  identify  software  and  system  design  considerations  for  an  extension  of  the  work 
in  a  Phase  II  effort. 

Results: 

1 .  A  general  review  of  digital  data  sources  was  completed.  This  information  was 
interpreted  to  provide  Phase  II  guidance. 

2.  Experimentation  with  selected  digital  data  sets  led  to  a  working  demonstration  of  the 
production  of  near  real  time  LOTS  planning  information  for  a  selected  area  near 
Camp  Le  Jeune,  NC.  The  data  along  with  operational  software  were  delivered  to 
the  SBIR  technical  monitor  at  USACEAVES. 

3.  "Processing  hypotheses"  were  tested  in  a  two-way  exchange  of  geographically 
encoded  information  between  a  geoprocessor  and  a  numerical  (surf  forecasting) 
model.  We  were  able  to  demonstrate  the  development  of  high  resolution 
bathymetry  from  multiple  input  sources  using  a  surface  building  application.  The 
bathymetry  was  then  processed,  in  what  we  consider  to  be  a  typical  MGIS- 
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numerical  model  interaction,  into  numerical  model  input.  The  numerical  (surf)  model 
was  run  and  its  output  was  imported  into  our  conceptual  MGIS  for  visualization. 


4.  We  identified  two  general  modes  of  MGIS  operation  as  the  basic  design  concept  for 
Phase  11 : 

a)  Climate  Mode,  in  which  knowledge  of  a  LOTS  area  is  developed  as  a  general 
review  of  trends  in  an  area’s  environmental  processes  (from  wind,  wave,  surf, 
current  and  tide  histories).  The  characteristic  data  for  this  mode  are  statistical 
reports  and  graphs; 

b)  Forecast  Mode,  in  which  observations  and  analysis  are  made  through 
visualization,  geoprocessing  and  numerical  model  interaction  in  a  detailed 
geospatial  context. 


Project  Narrative 

Our  View  of  GIS  Support  to  LOTS 

A  LOTS  planner  is  concerned  with  moving  material  and  people  from  sea  to  shore  using 
transports  and  assets  whose  capabilities  are  affected  by  environmental  conditions. 
Much  of  the  LOTS  planning  problem  focuses  on  gaging  vessel  and  asset  capabilities 
relative  to  conditions  predicted  for  many  types  of  environmental  conditions,  particularly 
waves,  surf,  currents,  tides  and  winds.  Other  types  of  information  such  as  water 
temperature,  air  temperature  and  bottom  characteristics  my  also  be  important. 
Bathymetry  is  of  key  interest  because  depth  variations  have  a  strong  affect  on  wave 
refraction  and  surf  characteristics. 

We  know  that  GIS  can  be  used  to  present  a  georeferenced  view  of  an  area  of 
operations  as  a  collection  of  thematic  maps.  A  thematic  map  can  represent  anything-  a 
detailed  shoreline,  bathymetric  contours,  predicted  wave  activity,  the  mapped  positions 
of  submerged  or  shore-based  structures,  etc.  In  our  current  GIS  design,  there  are  no 
effective  limitations  as  to  the  number  or  content  of  thematic  components.  The  only 
requirement  is  that  all  components  have  been  transformed  (geometrically  and 
geographically)  to  a  common  projection. 

A  series  of  standard  geoprocessing  functions  can  be  invoked  to  allow  an  operator  to 
"navigate"  the  GIS-  changing  viewing  opportunities  by  theme  or  scale,  calculating 
distances  and  bearings  from  point  to  point  or  initiating  queries  based  on  location  or 
(thematic)  attributes. 

At  a  general  level,  GIS  is  a  tool  that  can  already  provide  a  substantial  assist  to  LOTS 
planning  and  decision  support  by  visualizing  the  changing  features  of  the  nearshore 
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environment-  How  is  the  shoreline  oriented?  Where  are  significant  shorebased 
features?  What  are  the  observed  (mapped)  hydrologic  and  geologic  features?  We 
see  no  technical  problems  in  producing  a  visual,  geographic  data  base  of  this  class  of 
littoral  zone  data.  The  static  and  low  dynamic  thematic  components  of  a  MGIS  can  be 
developed  from  numerous  digital  and  analog  sources. 

The  question  is  how  to  augment  a  MGIS  with  processing  capabilities  from  the  predictive 
side  of  things-  where  maps  based  on  numerical  model  output  will  provide  a 
visualization  of  the  highly  dynamic  littoral  zone  at  times  and  scales  selected  by  the 
planner. 


MGIS  -  Numerical  Model  Interaction 

A  Phase  I  goal  was  to  explore  the  processes  through  which  static  and  low  dynamic  data 
could  be  located,  acquired  and  integrated  and  through  which  the  input  requirements  of 
a  numerical  forecast  model  could  be  developed.  We  selected  the  Navy  surf  model  as  a 
testbed. 

Model  Inputs 

Bathymetry 

Bathymetry  is  a  critical  element  in  many  forecast  models  which  MGIS  is  likely  to 
interact.  It  is  used  in  grid  cell  format  to  calculate  wave  refraction  in  near  shore  models 
(e.g.  STWAVE)  in  profile  format  for  surf  zone  calculations  and  surf  statistics. 

Sources  of  digital  bathymetric  data  were  determined  to  be  well-established  but  varying 
widely  in  geographic  coverage,  scale  and  accuracy.  From  the  civilian  perspective,  surf 
zone  bathymetry,  the  data  that  would  actually  be  required  for  surf  zone  forecasting,  is 
largely  non-existeiit,  particularly  in  areas  where  LOTS  exercises  may  occur.  Military 
beach  studies  are  the  main  source  of  very  high  resolution  surf  zone  bathymetry  at  a 
reasonable  number  of  locations. 

Expanding  GIS  to  include  tools  for  building  gridded  bathymetry. 

The  critical  importance  of  bathymetry  requires  that  MGIS  incorporate  tools  for  building 
bathymetric  maps.  In  Phase  I,  various  methods  of  creating  high  resolution  depth  charts 
from  digital  and  analog  sources  were  identified  and  investigated. 

The  basic  input  to  depth  chart  production  is  a  geoencoded  sounding-  a  longitude  (or 
projected  easting),  a  latitude  (or  projected  northing)  and  an  observation  of  depth.  Data 
in  this  form  is  available  in  several  digital  data  bases  (e.g.  NOS  Hydrographic  Data 
Base,  Navy  DBDB(n) ).  Data  can  also  be  produced  in  this  form  by  digitizing  nautical 
charts.  Given  the  convergence  of  DGPS  and  sensor  technology,  very  high  resolution 
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data  in  this  form  may  ultimately  be  developed  in  near  real  time  by  special  forces 
operations  before  amphibious  and  LOTS  operations. 

SURFACE  is  a  tool  developed  by  DDS  that  can  be  used  to  generate  gridded 
bathymetry  from  geoencoded  sounding  data.  Its  important  operational  aspects  are: 

1 .  Accepts  delimited  ASCII  input. 

2.  Automatically  transforms  input  projection  and  datum  to  a  user  defined 
"project"  projection  and  datum. 

3.  Operates  in  an  APPEND  mode  to  allow  multiple  source  inputs. 

4.  Permits  interactive  definition  of  spatial  limits. 

5.  Offers  a  selection  of  linear  or  non-linear  interpolation  schemes  (Trend 
Surface,  Nearest  Neighbor,  Inverse  Distance,  Exponential  Decay). 

6.  Produces  output  at  a  user  defined  spatial  resolution.  This  feature  is  very 
useful  in  relation  to  models  that  require  flexibility  in  defining  input  grid  sizes. 
Wave  refraction  models  are  an  example.  The  SURFACE  tool  can  be  used  to 
efficiently  build  variable  sized  output  grids  from  a  single  input  source. 

Figure  (1 )  is  an  example  of  SURFACE  output.  The  inputs  were  a  combination  of 
digitized  extractions  from  a  NOAA  chart  and  digital  STOIC  chart  soundings. 


figure  1 
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The  significance  of  this  analysis  is: 

1 .  The  product  is  a  thematic  component  of  a  MGIS  project. 

2.  The  product  is  a  geographically  accurate  visualization  of  depths  at  a  high 
resolution. 

3.  The  overall  accuracy  of  the  interpolated  surface  is  defensible  as  a  function  of  the 
accuracy  and  density  of  the  input  soundings.  The  interpolation  algorithm  effectively 
"honors  the  data." 

4.  The  product  is  in  a  format  that  is  compatible  with  near  shore  wave  refraction  models 
that  require  gridded  bathymetry  as  input. 

5.  The  source  of  the  data  is  not  critical.  Geoencoded  soundings  can  be  accessed 
from  existing  data  bases,  extracted  from  nautical  charts  or  collected  from  active 
sensors.  The  important  operational  factor  is  geography  and  data  format.  We  also 
note  that  bathymetry  can  be  progressively  refined  with  successive  data  inputs. 

6.  The  product  was  developed  in  a  few  minutes  of  processing  time  on  a  $3,000.00 
Pentium  PC. 

The  demonstration  of  the  applicability  of  SURFACE  to  MGIS-numerical  model 
interaction  is  considered  important  because  it  represents  a  capability  that  can  be 
leveraged  in  Phase  II  development  through  a  manageable  work  plan  effort. 

In  investigating  the  requirements  of  other  models,  it  was  noted  that  a  valuable  addition 
to  SURFACE  would  be  an  option  for  "data  smoothing".  Among  a  number  of 
approaches,  the  application  of  a  median  filter  (a  convolution  operation)  during  the 
surfacing  runstream  appears  to  be  an  effective  and  easily  implemented  solution. 

Our  investigations  showed  that  SURFACE  has  applications  potential  for  many  forms  of 
data  generation  in  MGIS-  numerical  model  interaction.  The  outputs  of  many  numerical 
models,  including  wave  refraction,  current  and  wind  forecast  and/or  hindcast,  are 
structurally  similar-  a  predicted  value  at  a  grid  location.  Any  data  structured  as  a  set  of 
coordinates  and  associated  values  can  be  processed  with  SURFACE. 

Tide  Adjustment 

GIS  tools  can  be  applied  to  make  a  "global"  adjustment  of  depths  in  the  surf  zone. 
Digital  addition  and  subtraction  can  be  operated  on  the  baseline  bathymetry  to  apply  a 
very  simple  correction  based  on  an  externally  derived  value  for  tidal  stage  at  a  fixed 
time. 


8 


Over  larger  areas,  modeling  tidal  effects  will  require  a  complex  MGIS-model  interaction 
in  order  to  visualize  the  variation  in  tide  levels  from  location  to  location. 

Shoreline  Shape  and  Orientation. 

The  LOTS  planner's  view  of  surf  zone  and  near  shore  bathymetry  must  be 
complemented  by  detailed  information  about  the  shape  and  orientation  of  the  shoreline. 
The  importance  of  shoreline  detail  is  that  it  provides  a  reference  for  the  definition  and 
orientation  of  surf  zone  bathymetry  profiles  and  wave  refraction  bathymetry  grids. 
Additionally,  shoreline  maps  will  typically  include  offshore  "blocking  features"  the  size 
and  position  of  which  can  limit  waves  from  reaching  specific  locations. 

Sources  of  georeferenced  shoreline  data  are  numerous  and  GIS  tools  for  processing 
shoreline  data  sets  are  well  developed.  From  the  civilian  perspective,  very  large  scale 
shoreline  coverage  is  severely  limited. 

In  MGIS,  "shoreline"  and  "detailed"  are  relative.  Shoreline  is  a  more  or  less  accurate 
representation  of  the  land/water  interface  (hopefully  compiled  relative  to  a  known  tide 
stage).  The  degree  of  detail  with  which  it  is  represented  is  a  function  of  the  scale  at 
which  it  was  mapped.  Figure  (2)  illustrates  relative  differences  in  a  representation  of 
the  same  shoreline  developed  from  two  sources  at  scales  of  1  ;2000000  (DMA's  Digital 
Chart  of  the  World)  and  1:100000  (USGS  DLG  III),  it  points  up  the  requirement  for 
detailed  coastal  maps  when  working  in  the  littoral  zone.  The  global  coverages  of  DCW, 
CIA  WDBII  are  developed  at  so  small  a  scale  as  to  be  unusable  in  the  very  nearshore. 


Initially,  the  bathymetric  map  (MGIS  SURFACE  output)  is  a  raster  image  that  contains 
no  detailed  shoreline  information.  A  vector  shoreline  can  be  laid  on  the  bathymetry  or, 
using  specialized  GIS  tools,  a  land  mask  can  be  created  to  overwrite  existing  depth 
values.  Figures  (3)  and  (4)  illustrate  each  approach. 
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figure  3  figure  4 


Once  bathymetry  and  shoreline  are  developed  as  thematic  components  of  our 
conceptual  MGIS,  the  basis  for  information  exchange  with  the  numerical  model  is 
established. 

Meteorological  and  Sea  State  Inputs 

In  addition  to  bathymetric  profiles,  the  SURF  model  requires  input  of  (deep  water)  wave 
height,  period,  direction  and  wind  speed  and  wind  direction  (preferably  from  local 
observation).  This  information  can  be  developed  from  hindcast,  forecast  or  near  real 
time  buoy  and  ship  observations  processed  through  regional  wind  wave  models 
(SHALWV  per  RTWF  document)  and  near  shore  wave  refraction  models  (STWAVE  per 
RTWF  document  and  HISWA  according  to  researchers  at  Woods  Hole). 

MGIS  could  rely  on  direct  links  to  Naval  regional  and  near  shore  wind  and  wave 
forecasts  (e.g.  NOGAPS,  NORAPS,  GOSWM,  WAM  area  all  evolving 
model/databases).  Because  model  output  is  geographically  and  geometrically 
consistent,  integration  and  visualization  of  gridded  wind  and  wave  field  data  is  implicit 
in  the  MGIS  design.  A  LOTS  planner  could  interactively  extract  forecast  wind  and  wave 
data  from  visualized  model  grids  and  pass  these  directly  to  wave  refraction  and  surf 
models. 


Feeding  the  SURF  Model 

Wind,  tide  and  sea  state  inputs  are  requested  from  the  LOTS  planner  interactively 
through  the  SURF  model  runstream.  These  are  individual  data  elements  which  we 
assume  the  LOTS  planner  will  have  preprocessed  and  on-hand  before  executing 
SURF.  Standard  GIS  inventory  and  query  tools  allow  very  rapid  development  of  this 


10 


class  of  inputs  from  available  georeferenced  data  or  from  external  models  (in  particular, 
Tide  level  forecasting). 

A  key  operational  issue  in  using  MGIS  to  feed  SURF  is  the  development  of  depth 
profile  input  from  GIS  bathymetry. 

The  depth  profile  required  by  SURF  is  a  collection  of  depths  and  distances  at  regular 
intervals  along  a  line  perpendicular  to  the  shore.  Two  methods  of  interactive  SURF 
profile  building  were  evaluated.  Both  assume  that  high  resolution  bathymetry  has  been 
(will  be)  developed  as  a  thematic  component  of  MGIS.  (We  note  that  the  SURF  model 
can  be  run  without  "real"  bathymetry  if  an  estimated  bottom  slope  is  supplied.  Using 
actual  bathymetry  is  preferred.) 


GIS  Profile 


We  have  developed  a  GIS  profiling  tool  that  extracts  data  from  a  GIS  theme  at  intervals 
along  a  user  defined  transect.  The  profile  extraction  can  be  exported  as  an  ASCII  file 
in  the  following  format: 

Easting  Northing  Longitude  Latitude  Depth 

286295.56  m.  3824740.34  m.  770  19’  43.65"  W  340  32'  38.29"  N  -0.0 

Positions  along  the  transect  are  reported  in  projected  and  geodetic  coordinates  from  an 
initial  to  terminal  point.  The  point  to  point  distances  required  by  the  current  version  of 
SURF  can  be  easily  developed  from  the  basic  profile  data. 

The  profile  is  developed  visually  by  point  and  click  on  GIS  bathymetry  and  shoreline.  A 
number  of  profiles  can  be  generated  in  rapid  succession.  The  profile  interface  provides 
the  user  with  information  on  transect  length  and  azimuth  so  the  user  has  latitude  in 
defining  distances  and  bearings  relative  to  a  visual  presentation  of  shoreline. 


11 


1  GIS  Distance 

Projection  : 

UNIVERSAL  TRANSVERSE 

Northing: 

3823108.23m 

Easting: 

286874.71m 

Latitude: 

34*  31-  45.78"  N 

Longitude: 

77*  19- 19.48“  W 

Pixel  Value  : 

-3.9988 

Initial  Lat: 

34*  32‘  04.38"  N 

Initial  Lon: 

77*  19*  36.80“  W 

Forward  Azimuth:  1 42.3760^ 

Backward  Azimuth:  322.3787* 

Distance: 

0.3908nmi 

T otal  Distance: 

0.3908nmi 

Left  mouse  to  select  position 
Right  irtouse  to  exit 

Value  :  -5.108497620 

Easting;  287208.02  m  Latitude:  34*  31*  39  85"  N 
{  Northing:  3822917.76  m.  Longitude:  77*  1 9' 06.24"  W 
I  Exit  I  Report  j  Help  j 


Operationally,  the  user  has  the  option  to  build,  in  a  more  or  less  free  form  session,  a 
variably  spaced  grid  based  on  his  identification  of  areas  of  interest  in  the  surf  zone.  A 
collection  of  profiles  can  be  successively  processed  by  SURF. 

CIS  Grid  Sampling  of  Bathymetry 

It  was  observed  that  however  useful  interactive  profiling  may  prove  to  be,  a  requirement 
exists  to  develop  profiling  by  laying  a  rotated  grid  on  bathymetry  and  extracting  depths 
automatically  at  grid  intersect  points.  We  have  developed  a  GIS  tool  for  over  the  past 
fourteen  years  for  customers  requiring  grid  sampling  that  can  be  adapted  to  meet  this 
requirement.  Our  testing  was  limited  to  a  concept  demonstration  of  grid  definition  and 
grid  rotation. 

The  SAMPLING  tool  is  operated  on  a  visual  display  of  bathymetry  and  shoreline. 
Horizontal  and  vertical  grid  spacings  and  grid  orientation  are  defined  interactively  by 
the  user.  The  grid  overlay  is  visual  and  can  be  adjusted  as  necessary. 


We  believe  that  this  sampling  tool  can  be  developed  as  the  primary  method  of 
preparing  bathymetry  profile  and  bathymetry  grid  input  to  various  numerical  models  by 
average  users.  Its  principal  advantage  is  the  ease  and  consistency  with  which  an 
"extraction  grid"  can  be  developed. 

In  Phase  II.  with  variable  resolution  grids  of  wind  and  wave  field,  bathymetry  and 
shoreline,  the  grid  sampler  should  be  able  to  extract  all  the  local  parameters  required 
by  SURF.  An  exception  would  be  tide  heights.  These  must  be  developed  externally 
through  a  tide  forecast  model. 

The  expansion  of  the  “extraction  grid”  capability  would  involve  building  a  capability  to 
extract  multiple  parameters  (from  a  layering  of  thematic  components)  at  each  grid 
center  or  a  user  defined  grid  corner. 


Visualizing  SURF  Model  Output 

The  PROFILE  tool  was  used  to  build  four  depth  profiles  in  the  surf  zone.  The  initial 
ASCII  output  was  processed  into  a  listing  of  distances  and  depths  as  required  by 
SURF.  The  files  were  transferred  via  e-mail  to  NSI  offices  where  the  SURF  model  was 
run.  SURF  output  was  returned  via  e-mail.  We  believe  that  this  transaction  was 
representative  of  a  real  world  situation  in  which  two  LOTS  operators  are  working 
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together  via  a  digital  comm-link  (clearly,  faster  and  more  secure  communications  than 
e-mail  vi/ould  be  used). 


SURF  output  includes:  the  distance  from  shore  of  a  modeled  observation  and  the 
water  depth,  significant  breaker  height,  maximum  breaker  height,  percent  breaking 
waves,  wave  length  and  littoral  current  at  that  distance. 

Geoprocessing  of  SURF  output  first  requires  a  conversion  of  distances  from  shore  to 
coordinate  values.  In  this  test,  the  conversion  was  accomplished  with  Excel  using  its 
trig  functions.  X  and  Y  are  the  computed  UTM  Zone  18N  coordinates  (NAD  27). 


Index 

Distance  Depth 

Bht 

Max 

%B 

Wave  L 

Current 

X 

Y 

21 

6982.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.08 

288711.3 

3822690 

22 

6972.3 

17.2 

3.1 

4.8 

0.1 

192.2 

p.08 

288709.4 

3822692 

23 

6962.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288707.6 

3822695 

24 

6952.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288705.8 

3822697 

25 

6942.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288703.9 

3822700 

26 

6932.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288702.1 

3822702 

27 

6922.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288700.3 

3822705 

28 

6912.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288698.4 

3822707 

29 

6902.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288696.6 

3822710 

30 

6892.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288694.8 

3822712 

31 

6882.3 

17.2 

3.1 

4.8 

0.1 

192.2 

0.09 

288692.9 

3822714 

32 

6872.3 

17.2 

3.1 

4.7 

0.1 

192.2 

0.09 

288691.1 

3822717 

33 

6862.3 

17.2 

3.1 

4.7 

0.1 

192.2 

0.09 

288689.3 

3822719 

34 

6852.3 

17.2 

3.1 

4.7 

0.1 

192.2 

0.09 

288687.4 

3822722 

35 

6842.3 

17.1 

3.1 

4.7 

0.1 

191.6 

0.09 

288685.6 

3822724 

36 

6832.3 

16.9 

3.1 

4.7 

0.1 

191 

0.09 

288683.8 

3822727 

37 

6822.3 

16.8 

3.1 

4.7 

0.1 

190.5 

0.09 

288681 .9 

3822729 

38 

6812.3 

16.8 

3.1 

4.7 

0.1 

190.5 

0.09 

288680.1 

3822731 

39 

6802.3 

16.8 

3.1 

4.7 

0.1 

190.5 

0.09 

288678.3 

3822734 
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The  requirement  to  “build  in”  georeferencing  points  to  an  important  problem  in  MGIS- 
model  interaction.  Surf  model  output  is  not  explicitly  geographic.  In  the  absence  of  a 
mechanism  for  calculating  and  appending  geographic  coordinates  to  each  record  of  the 
surf  model  output,  the  results  (in  column  2)  are  simply  point-to-point  distances  in  an 
abstract,  “model  space”.  In  MGIS-model  interaction,  “model  space”  and  geographic 
space  must  converge  before  model  output  can  be  visualized  in  a  geospatial  context. 

Because  the  profile  (model  bathymetry  input)  is  extracted  from  a  georeferenced  source, 
we  were  able  to  recover  geographic  coordinate  values  for  each  surf  model  element 
based  on  profile  start  and  stop  locations  and  transect  angle. 

If  the  surf  model  is  to  be  made  available  to  MGIS  operators,  this  convergence  problem 
will  require  a  transparent  solution. 

SURF  model  output  in  this  final  form  (*.csv)  is  ready  input  to  the  SURFACE 
geoprocessing  function.  The  initial  plot  of  surf  model  output  is  illustrated  below  in  the 
following  figure.  The  “lines”  are  in  fact  points  representing  the  individual  model 
observations  at  calculated  distances  from  shore.  There  were  four  input  profiles  to  the 
model  and  four  output  sets. 
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Applying  the  SURFACE  process  to  SURF  output  results  in  the  following  visualization  of 
significant  wave  heights  in  the  surf  zone  (scale  is  1 :17000)  in  feet: 


Cotoi  Bar 


surfmodS  Field  5 


m. 


0.28  0.87  1.46  2.05  2.64  3.23  3.82  4.41  5.00 


The  preceding  image  is  important  from  several  aspects; 

1 .  In  our  view,  it  objectively  demonstrates  the  feasibility  of  integrating  geoprocessing 
with  numerical  modeling  in  support  of  LOTS.  Model  input  bathymetry  was  developed 
inside  our  conceptual  MGIS  from  a  STOIC  chart  base.  Model  input  was  developed 
using  an  easily  operated  GIS  tool.  Model  output  was  easily  and  quickly  ingested  and 
processed  in  a  single  operation  (SURFACE). 

The  result  is  a  georeferenced  view  of  predicted  wave  activity  in  the  surf  zone  near  New 
River  Inlet,  NC.  It  functions  as  a  highly  specialized  nautical  chart. 


2.  The  image  is  a  map-  a  scaled,  projected  representation  of  features  in  the  surf  zone. 
Using  it,  a  LOTS  planner  can  begin  to  determine  what  areas  of  the  surf  zone  landscape 
are  predicted  to  be  safe  for  transit  by  specifically  configured  vessels. 

3.  The  image  is  a  georeferenced,  thematic  component  of  a  concept  MGIS.  Its  spatial 
relationship  (coincidence,  proximity,  adjacency)  with  other  MGIS  themes  is  implicit. 

4.  The  process  can  be  repeated  through  different  times  to  develop  a  series  of  surf 
zone  forecasts  as  MGIS  components.  The  MGIS  can  have  a  temporal  as  well  as  a 
spatial  and  thematic  dimension.  Basic  GIS  theory  implies  an  ability  to  build  a  map  of 
composite  wave  activity  through  several  forecast  periods  for  example. 

Further  MGIS  Analysis 

The  model  output  is  a  thematic  component  of  our  conceptual'  MGIS-  the  theme  in  this 
case  is  significant  wave  height.  Additional  analysis  can  be  performed  on  the  data  using 
MGIS  tools.  For  example,  basic  statistics  on  the  frequency  and  spatial  distribution  of 
significant  wave  heights  can  be  developed  by  processing  the  thematic  map  through  a 
statistical  summary  and  display  tool.  The  following  figure  illustrates  the  output  of  this 
tool  as  it  was  run  on  the  model  output: 


The  left  panel  illustrates  a  graph  of  frequency  distribution  for  various  wave  heights 
within  the  model  viewed  in  the  right  panel.  The  statistics  area  shows  the  value  of  the 
current  highlighted  area  (white  contour)  as  3.03  feet.  The  mean  wave  height  is  given 
as  2.06  feet  and  standard  deviation  within  the  model  area  is  0.9194.  As  the  operator 
steps  through  the  graph,  the  image  is  instantly  updated  to  display  the  spatial 
distribution  of  the  current  step  value. 

MGIS  operations  such  as  this  (Object  Info)  have  value  as  a  support  to  basic 
visualization.  We  anticipate  an  adaptation  and  expansion  of  these  capabilities  in 
Phase  II. 


Data  Issues  and  Observations 


General  Listing  of  LOTS-critical  parameters 

We  estimate  that  there  are  as  many  as  fifty  environmental  parameters  with  significance 
to  LOTS  that  must  be  processed  in  a  geospatial  context  either  as  thematic  maps  in  a 
forecast  mode  or  as  mixed  media  objects  in  a  climate  mode. 

The  following  list  outlines  those  identified  in  Phase  I.  Their  major  sources  are  detailed 
in  an  appendix  to  this  report. 

Climatic  parameters 

Wind  (speed,  direction,  time)  and  (percent  occurrence  for  speed  vs.  direction,  time  period) 
Precipitation  (type,  amount,  time  period)  and  (percent  occurrence) 

Severe  storms  (category,  track,  time,  percent  occurrence) 

Visibility  (distance,  time  of  day,  percent  occurrence,  time  period) 

Cloud  cover  (amount,  time,  probability  of  occurrence) 

Barometric  pressure  (pressure,  time) 

Moonlight  (phase,  time,  position) 

Sunlight  (sunrise,  sunset,  twilight) 

Superstructure  icing  (possibility  of  occurrence) 

Wind  chill  (temperature,  time) 

Oceanographic  parameters 

Waves  (height,  period,  time)  and  (probabilities  of  spectrum) 

Currents  (type,  speed,  direction,  time,  time  period) 

Astronomic  tide  (type,  ranges,  height,  time,  coefficients,  sample  curve) 

Non-tidal  fluctuations  (height,  time) 

Sea  ice  (thickness,  percent  cover,  time,  time  period) 

Water  temperature  (  percent  occurrence) 

Physical,  geological,  biological  features 
Shoreline 

Bathymetry/topography 

Beach  profile  (depths  across  surf  zone  onto  beach) 

Surface  materials  type  (thickness,  type) 

Areas  of  high  sedimentation  or  erosion  (area,  type) 
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Landform/ecology/land  use  (type) 

Vegetation  (type) 

Rivers,  lakes 
Watershed  gradients 

Other  Physical  and  Cultural  Features 
Buildings 
Roads 
Piers,  docks 
Navigation  aids,  buoys 
River  mile  points 
Airstrips 

Misc.  man-made  structures 

Landmarks 

Political  boundaries 

Place  names 

Vessel  traffic 

Obstructions 

Operational  parameters 
Commercial  shipping  lanes 
Navy  pre-planned  "Q"-routes 
Landing  approaches 
Penetration  points  and  lanes 
Target  areas,  command  locations 
Passage  height  under  bridges 
Landing  craft  operational  parameters 
Site  objectives 

High  Resolution  Bathymetry  and  Shoreline 

Although  Work  in  Phase  I  was  undertaken  from  a  completely  civilian  perspective,  it  was 
supported  by  reasonably  good  data.  We  were  able  to  locate  accurate  unclassified 
bathymetry  (STOIC)  and  shoreline  data  in  digital  form.  Sea  state  and  meteorological 
observations  were  available  in  near  real  time  from  a  number  of  weather  and  sea  state 
observation  facilities  based  on  and  off-shore  (OMAN  and  NDBC  stations). 

LOTS  planning  operations  will  most  likely  not  have  immediate  access  to  high  resolution 
digital  data,  therefore,  a  requirement  for  CIS  tools  will  be  developed  that  focus  on 
independent  data  creation  by  a  LOTS  planner. 

GIS  tools  currently  exist  which  are  capable  of  supporting  high  resolution  data 
development  in  key  areas  such  as  bathymetry  and  shoreline  mapping.  We  have 
developed  and  demonstrated  data  generation  capabilities  which  can  be  integrated  in 
MGIS: 

1 .  Table  digitizing  and  heads  up  screen  digitizing  which  will  allow  a  LOTS  planner  to 
generate  geoencoded  soundings  from  digital  nautical  charts.  (Note  that  creating  a 
georeferenced  digital  nautical  chart  can  be  accomplished  by  standard  scanning  and 
geoencoding  methods  in  a  few  hours.  This  is  a  mature  technology.) 
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2.  Combining  remote  sensing  and  heads  up  screen  digitizing  which  allow  a  LOTS 
planner  to  generate  high  resolution  shoreline  maps  from  geoencoded  satellite  or 
airborne  imagery.  We  believe  that  it  is  critical  to  recognize  that  geoencoded 
imagery  is  a  map  that  can  be  used  in  MGIS  as  the  basic  visualization  context  for 
LOTS  planning  operations. 

3.  GPS-based  data  logging  which  can  be  used  in  field  studies  for  the  development  of 
point  observations  of  virtually  any  measurable  event  (e.g.  water  temperatures, 
clarity,  salinity,  bottom  materials,  plant  and  animal  communities,  etc.).  We  believe 
that  GPS-based  data  logging  will  become  a  standard  tool  for  high  resolution  field 
reconnaissance  and  mapping. 

Regional  and  Global  Data 

LOTS  planners  will  often  rely  on  medium  and  low  resolution  data  for  generalized 
mapping  and  interpretation  support. 

in  our  Phase  I  investigations  of  regional  and  global  data,  we  surveyed  the  developing 
information  systems  at  Fleet  Numerical  and  Navoceano  (e.g.  Navoceano's  IDBMS). 

The  data  holdings  and  basic  data  characteristics  are  well  documented  in  literature 
distributed  by  COMNAVMETOCCOM. 

Phase  I  has  demonstrated  the  general  ability  to  process  and  to  visualize  geoencoded 
data.  Most  of  the  data  referenced  in  the  Meteorological  and  Oceanographic  Data  for 
Modeling  and  Simulation  publication  is  either  geoencoded  point  or  gridded  data.  If  we 
can  extract  data  subsets,  we  can  process  them  into  a  geographic  visualizations  (maps) 
and  integrate  them  as  thematic  components  of  a  MGIS  project.  Indications  from  system 
descriptions  of  Navoceano  IDBMS,  NMDC  GLIS  and  NOAA's  NODDS  are  that  data 
subsets  can  be  extracted.  The  same  geoprocessing  tools  that  were  applied  to  high 
resolution  data  can  be  used  with  data  developed  from  regional  and  global  data.  From 
the  standpoint  of  projected  MGIS  capabilities,  the  only  operational  difference  between 
the  geoencoded  STOIC  soundings  and  data  from  DBDB-1,  DBDB-2,  DBDB-5  or 
ETOPO  is  scale. 

SQL-based  mechanisms  for  data  set  extraction  are  being  developed  within 
Navoceano's  IDBMS.  At  this  point,  we  believe  that  an  efficient  course  for  Phase  II  work 
would  be  to  use  IDBMS  as  a  gateway  for  regional  and  global  data  retrieval  in  the 
development  of  MGIS  projects. 
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Special  Issue:  The  Presentation  of  Tabular  and  Unmapped  Imagery 

Not  all  LOTS  operators  will  be  directly  concerned  with  mapping  and  numerical  model 
interaction.  There  is  a  class  of  operators  who  are  concerned  with  rapidly  developing 
basic  site  descriptions  in  tabular  and  graph  form. 

Geoprocessing  technology  has  advanced  to  a  point  at  which  geoencoded  objects  can 
be  attributed  with  external  objects  or  applications.  It  is  currently  possible  to  link  an 
object  or  location  in  a  map  with  a  document,  spreadsheet,  static  image  or  mpeg. 

The  possibilities  for  extending  MGIS  to  operation  in  a  historic  or  "climate  mode"  are 
very  interesting.  We  believe  that  a  basic  map  interface  can  be  developed  as  a  set  of 
visual  access  keys  to  tabular  and  non-georeferenced  imagery.  A  user  can  click  on  a 
feature  in  a  map  and  immediately  access  an  mpeg  player  to  run  video  of  an  aerial 
overflight  or  to  view  a  graph  of  historic  wave  activity  or  climate  data. 

The  conceptual  MGIS  demonstrated  this  operating  concept  by  accessing  an  external 
tide  prediction  tool  within  the  general  MGIS  application  when  the  user  clicked  anywhere 
in  the  offshore  area  of  one  of  the  thematic  components. 

Without  trivializing  the  data  development  and  data  management  challenges,  we  note 
that  "mixed-media"  capabilities  are  available  in  current  GIS  design.  We  believe  that 
they  can  perform  valuable  functions  within  MGIS. 


Summary  and  Conclusions 

1 .  Geoprocessing  (GIS)  is  a  mature  technology  that  can  effectively  support  LOTS. 

This  Phase  I  effort  has  effectively  demonstrated  the  capability  of  PC-based 
geoprocessing  to  integrate  environmental  data,  including  the  output  of  a  numerical  surf 
forecast  model,  in  a  visual,  georeferenced  information  system.  The  general  feasibility 
of  MGIS  was  established  in  Phase  I. 

2.  A  critical  question  was  the  feasibility  and  degree  of  difficulty  involved  in  developing 
an  interaction  between  numerical  models  and  GIS.  To  the  extent  that  we  succeeded  in 
designing  a  basic  interaction  between  the  Navy  surf  model  and  a  conceptual  MGIS, 
general  feasibility  was  addressed. 

We  believe  that  the  MGIS-surf  model  interaction  demonstrated  a  basic  transaction 
structure  that  will  be  repeated  in  interactions  with  other  wave  and  surf  forecast  models, 
notably  the  current  and  future  generations  of  STWAVE  and  HISWA,  We  expect  to  offer 
LOTS  planners  a  range  of  model  interaction  options  based  on  their  level  of  comfort  with 
a  given  model. 
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We  are  convinced  that  the  “user  intensity”  of  MGIS/numerical  model  interactions  can 
be  managed  through  the  development  of  application  specific  software  in  a  Phase  II 
effort. 

3.  Phase  I  also  demonstrated  the  facility  with  which  critical  high  resolution  bathymetry 
can  be  developed  using  standard  geoprocessing  tools.  Several  methods  for 
developing  high  resolution  bathymetric  and  shoreline  maps  were  identified.  We  note 
that  our  experience  with  these  methods  is  significant  and  predates  this  Phase  I  work. 

4.  The  mechanics  of  populating  a  MGIS  project  with  geospatial  data  are  well 
established  in  current  GIS  technology.  Questions  concerning  actual  data  availability, 
quality  and  utility  in  LOTS  planning  remain.  We  see  major  issues  relating  to  the 
resolution  of  regional  and  global  data  sets  and  their  ultimate  utility  to  LOTS  planning  in 
the  littoral  zone. 

Some  50  LOTS-critical  data  requirements  were  identified  and  sources  of  these  data 
were  surveyed. 

The  majority  of  these  data  exist  as  geoencoded  grids  or  point  observations.  As  such, 
they  are  available  to  GIS  integration.  The  general  operating  theory  of  information 
systems  such  as  NOAA  NODDS  and  NAVOCEANO  IDBMS  indicates  that  global  and 
regional  data  location  and  retrieval  are  feasible. 

Because  much  of  the  data  is  distributed  in  ASCII  format,  data  ingest  issues  will  be 
manageable.  Data  in  this  form  can  be  processed  using  currently  designed 
geoprocessing  tools.  The  SURFACE  application  demonstrated  a  well  developed 
capability  for  ingest  and  processing  of  this  significant  class  of  geospatial  data. 
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